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Abstract 

We analyze here a model for single-electron charging in semiconductor quan- 
tum dots that includes the standard Anderson on-site repulsion (U) as well 
as the spin-exchange (Jd) that is inherently present among the electrons oc- 
cupying the various quantum levels of the dot. We show explicitly that for 
ferromagnetic coupling (J^ > 0), an s-d exchange for an S=l Kondo problem 
is recovered. In contrast, for the antiferromagnetic case, < 0, we find that 
the Kondo effect is present only if there are an odd number of electrons on 
the dot. In addition, we find that spin-exchange produces a second period in 
the conductance that is consistent with experimental measurements. 
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When a gate voltage is applied to a nano-scale semiconductor inversion layer (or quan- 
tum dot), electrons will flow one at a time across this device provided that the applied 
voltage is an integral multiple of the capacitance charging energy of the quantum dot. Ex- 
periments illustrating the principle of charge quantization by virtue of the charging energy 
have been performed recently on numerous semiconductor [IJ |2| [|3| j3j |J as well as super- 
conducting H nano-structures. We focus here solely on the semiconductor devices. It is 
now well-accepted || that in semiconductor quantum dots, the dominant contribution to 
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the capacitance charging energy, E c = ^ arises from the on-site Coulomb repulsion. Here 
C is the capacitance between the quantum dot, the tunnel junctions, and the electrical leads 
connected to the dot. Transport in quantum dots will be Coulomb limited if ksT < Ec and 
UbT > Ae, where Ae is the spacing between the single particle states of the dot. 

Because of the central role played by on-site Coulomb repulsions in the transport prop- 
erties of quantum dots, it is natural to model a quantum dot with a Hubbard-like model. 
In so far as a quantum dot can be reduced to a single site |7J with a charging energy U, the 
Anderson model || for the interaction of a magnetic defect coupled to a non-interacting sea 
of conduction electrons is appropriate [Tj: 

Ha = Y1 e k a L a ka + € d a da a da + ^ V kdWka a da + a da a ka) + U n d ^n di (1) 
k,a o" k,a 

= H + Vkd(al a a da + a da a ka ) (2) 

k,a 

where e d is the defect energy of the magnetic impurity, \%d the overlap integral between a 
band state with momentum k and the impurity, a k creates an electron in the band states, a) d(T 
creates an electron with spin a on the impurity, and n da = a\ a a dcT is the number operator 
for an electron of spin a. As a consequence of the on-site repulsion, the single particle 
states on the impurity have energies, e d and e d + U. At high temperatures, the density of 
states of this model has two Lorentzian peaks centered at these two energy levels. At low 
temperatures, however, the Anderson model displays a Kondo resonance JT0[] at the Fermi 



level. Although the Kondo resonance is expected to occur for any value of the defect energy 
within the range —U < e d < 0, it is most favourable at the defect energy corresponding 



to the greatest stability of the local moment at the d- impurity, namely, = — At this 
energy Ha is particle-hole symmetric, and the Kondo resonance is pinned at ef = 0. The 
single-particle states e d and e d + U lie symetrically, then, around the Fermi energy. When 
the chemical potential of the source lead coincides with the energy of the Kondo resonance, 
a single electron should charge the dot. In the symmetric limit, this state of affairs should 
obtain at half-integer multiples of the charging energy. Thus far, no experimental hint of 
the Kondo resonance has been observed in quantum dots in zero bias voltage. 

It is precisely the conditions under which the Kondo effect should be observed in quantum 
dots that we address here. While it may be premature to draw any conclusion from the lack 
of experimental confirmation of the Kondo effect, it is certainly appropriate to investigate 
the validity of the Anderson model to a quantum dot. It is in addressing this issue that we 
are 1) able to predict even-odd charging effects in semiconductor quantum dots as well as 
2) a suppression of the Kondo effect when a quantum dot has an even number of electrons. 
The most obvious inadequacy of the Anderson model in the context of quantum dots is the 
truncation of the multiple electronic levels on a quantum dot to a single state. If multiple 
electronic levels are included on the dot, then other Coulomb interactions besides the on-site 
U become relevant. A key quantity that comes into play is the intrinsic Coulomb exchange 
energy, J d between two levels. Consider for the moment a two-orbital model for a quantum 
dot. For two degenerate levels and U ~ the direct Coulomb exchange integral, the energy 
of the two-body states predicted by this model are 2e d + U — ^f, 2e d + U, and 2e d + U + ^-r L . 
The energy of the 3-electron state is 3(6^ + U). Consequently, the charging energy depends 
on J d . In fact, the general role of spin-exchange is to introduce a spin-dependent charging 
energy that is determined by the parity of the total number of electrons on the dot. It 
is worth pointing out that charging experiments on controlled-barrier atoms in strong 
magnetic fields display systematic oscillation in the peak heights, widths as well as in the 
separations that are consistent with a second period in the conductance as a function of the 
applied gate voltage. Such systematic deviations have been attributed to a splitting of the 
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energy between the up and down Landau levels, rather than to spin exchange. We propose 
here that such trends are also consistent with a spin-exchange model. 

We explore then the simplest model of a quantum dot that includes the effects of spin- 
exchange. A natural way of including spin exchange is simply to introduce another level 
into the Anderson model. The only qualitative change this level is going to provide is 
the spin-interaction with the d-level. Consequently, we treat this level as a local spin, S. 
Because spin-exchange plays no role if S — 0, we will consider only the case in which the 
S-level is singly occupied, or equivalently, S z = ±|. Hence, large N expansion techniques 
are inappropriate to solving this problem. If we label the spin on the d-level of the impurity 
with Sd, we find that our Hamiltonian can be written as 

H = H A + Hj (3) 
Hj = -J d S d • S (4) 

s d = ^Kt - n di), S£ = 4 T a di , = a) di a dV (5) 

The first question we answer with this model is, does the Kondo effect still occur. Be- 
fore rigorous calculations are performed, a heuristic answer can be put forth immediately. 
Without loss of generality, the defect energy can be taken to be \e d \ U >> \J d \, as J d 
is typically a fraction of U. The tunneling rate to the dot is determined by the matrix el- 
ement Vk d - Because this quantity is an adjustable parameter determined by the width of 
the tunnel junction connecting the dot to the source lead, we can set \ J d \ > \Vk Fd \. This 
inequality is crucial in the analysis of what follows. For example, in the absence of the J d 
coupling, the form of the antiferromagnetic interaction that gives rise to the Kondo effect 
scales as \Jk\ ~ ^ Vk ^ . However, in the limit that | J d \ > \V kFd \, \ J d \ » \Jr\- That is, the 
exchange coupling exceeds the Kondo coupling and could hence ultimately conspire to mask 
the Kondo effect. Consider the case in which the d-level and the S-level are singly occupied. 
In the ferromagnetic case, J d > 0, the ground state of the dot is a triplet. A Kondo effect 
should result in this case that is determined by the total spin on the dot. However, in the 
antiferromagnetic case, J d < 0, the ground state on the dot is a singlet. As a consequence, 
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there is no net spin to couple to the conduction electrons and the Kondo effect is suppressed. 
In the antiferromagnetic case, there must be an odd number of electrons on the dot for the 
Kondo effect to be observed. This is the essential physics of this model. 

To prove the heuristic arguments given above, we diagonalize H in the subspace of 
all singly-occupied states on the dot. We first note that because [H ,Hj] = 0, we can 
work entirely with the eigenstates of the dot. Let us define a generalized eigenstate \q) = 
\N: Stot, <Stot)> where N refers to the number of electrons on the d-level of the dot. Recall, 
we have set S z = ±|. There are 8 eigenstates in the dot basis: |0; \, S^ ot = ±|) with total 
energy E=0 , the singlet |1; 0, 0) with energy e s = e d + ^p, the triplets |1; 1, S^ ot = 0, ±1), 
with energy e t = e d — and the doubly-occupied state, |2; |, S^ ot = ±|) with energy 2e d + U. 
Each of the Fermion operators as well as the bi-linears such as S d • S can be expressed in 
terms of these 8 basis states. Once this is done, matrix elements among these states can 
be determined straightforwardly. For an arbitrary wavefunction we are interested in 
solving the Schroedinger equation 

J2(q\H\p){ P m = E{q\i;) (6) 
p 

in the subspace of singly occupied states, N=l. This is the relevant phase space for consid- 
ering the Kondo effect. To reduce the full 8 x 8 to a 4 x 4, we rewrite the matrix elements 
involving the empty and doubly-occupied states in terms of the singly occupied states. The 
exact result is a 4 x 4 Hamiltonian matrix 



H= 







R n 




Ht,o 


- R n 






Ht,i 


- R n 


- R n 






HT 



V 

where the singlet and triplet Hamiltonians are 



t-i J 



(1;0,0|V> 

(i;i,o|V> 
<i;i,i|^> 

<1; 1, — 1|V> 



(7) 



H S = H C + e s + J2Qit\ H t , = H c + e t + J2Q 



(+) 



Ht,i = H c + e t + 2Q^\ H t -\ = H c + e t + 2<5^ 



(8) 
(9) 
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and Q and R are the matrix elements 

<& = ^ E V k >V k *[±(E + e k , -H c - 2e d - Ur l a + kla a ka - (E - e k - H^a+^a^] (10) 
1 kk' 

R *S = 4? E V*Vkl(E + *v -H e - 2e d - U)- 1 + (E - e k - H^a^a^ (11) 
v2 kk , 



The matrix elements R and Q contain all powers of the coupling to the conduction 
electrons. To lowest order, they scale roughly as ^fj-- The Hamiltonian matrix can be 
partitioned into lxl singlet and 3x3 triplet subspaces provided that the differences 
between the diagonal elements exceeds the off-diagonal matrix elements R and Q. The 
diagonal elements differ by the spin-exchange J d - Consequently, the partitioning into singlet 
and triplet subspaces is valid provided that | J d \ > The effective Hamiltonian in each 

subspace that is valid to second order in the coupling to the leads can be obtained by setting 
E = e s + H c and E = e t + H c in the denominators of R and Q and transcribing the basis state 
representation back to the original Fermion operators. In the ferromagnetic case (J d > 0), 
the reduced Hamiltonian in the triplet subspace is 

H trijlet = Hc + et+ J2 WlM^k) ~ E 4k'(4'^k) ■ Stot (12) 
kk'a kk' 

where H c is the Hamiltonian for the free conduction electrons, a is the Pauli spin matrix, 
ip k is the two-component spinor 



i>k 



the antiferromagnetic coupling constant is 

4k' = W k , d V k \ (- ^ 1 TT ■ - - - J— ) (13) 



e k > - 2e d - U + e t e k - e t 



and 



Wt* = \v»vJ \ u+ +-L-). (M) 



That Jl k , is negative can be seen immediately because the largest energy scale in the de- 
nominator is U which enters with a — sign. The spin interaction obtained in this limit is 
identical to the usual Kondo coupling except in this case the total spin on the dot enters. 
For a triplet state S to t = 1- Consequently, ferromagnetic exchange gives rise to a S=l Kondo 
problem. The S=l Kondo problem is an example of an undercompensated spin problem in 
which the conduction electrons only partially screen the spins on the dot. The remaining 
unscreened spin couples ferromagnetically to the conduction band. The only qualitative 
difference between the S=l Kondo problem and S=| is the behavior of the magnetic sus- 
ceptibility. As a result of the undercompensation, the susceptibility does not vanish at T=0 
in the S=l problem. 



Consider now the more experimentally-relevant antiferromagnetic case ||11|| . We recover 
in this limit an effective Hamiltonian of the form 

H spiet = Hc + €s + Yw^Ut) + 4 lK Jk 2k 'M[ a M ■ (4^) (is) 

fcfc' d k lk [ 

where W£ k , and J| fc/ are identical to their triplet counterparts with e s replaced by e t . As is 
evident the spin coupling only involves the conduction electrons and is 0{{Vkd) i )- Further, 
the overall sign of this interaction is negative or ferromagnetic. Consequently, there is no 
antiferromagnetic exchange interaction that can produce a Kondo effect to fourth order in 
the coupling to the band electrons. The physical origin of the absence of the Kondo effect 
here is the stability of the singlet on an energy scale J^. As a result, the Kondo coupling 
constant must be cut off at this energy scale. Consequently, it cannot diverge and give 
rise to a bound state at the dot. This result is consistent with the heuristic arguments of 



Ng and Lee on the role of spin exchange in a quantum dot and a mean-field limit of 



the 2- impurity Anderson model in the presence of spin-exchange ||13|| . Jones, Kotliar and 
Millis |TJj found that in the iV— > oo limit of this model, a phase transition occured which 
suppressed the Kondo effect if the bare exchange interaction exceeded a critical value. The 
critical condition is similar to the one used here, namely \J$\ > . There still remains 
one chance for the Kondo effect to be observed when < 0. If the number of electrons on 



the dot is odd, or equivalently we restrict ourselves to the N=0 subspace, the standard S=| 
Kondo problem is recovered. If experiments are going to detect the Kondo effect, the total 
number of electrons on the dot must be carefully controlled. 

In deriving the results in the Kondo regime, we have performed 2nd-order perturbation 
theory in the coupling to the leads. It is possible to construct a Schriefer- Wolff- type || 
transformation that eliminates the coupling to the leads. The result 



S 



2S + 



k,cr,a' 



'S + l 



s 



+2(S-s 



E, 



El 



(-) 



(-) 



E 



E, 



(-) 



(-) 



v (S + l S 
l-n d - a )+ I —r?r + 



(1 - n d - a , 



.E. 



(+) 



E\ 



(+) 



Ud-a' 



clod*' 

-h.c. (16) 



can be used to derive the effective Hamiltonians in the singlet and triplet subspaces. In the 
above, Eg t ' = — e s>t and E^ = — 2e d — U + e Sjt - In the singlet subspace, S = and only 
those virtual transitions involving the triplet state survive. This transformation successfully 
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eliminates the coupling to the leads in the limit « 1 and hence is consistent with the 
perturbative treatment developed here. 

We now calculate the conductance at finite temperature in the presence of spin-exchange. 
To facilitate this we need the average occupancy on the dot (n da ). This quantity is obtained 
by integrating the imaginary part of the d-electron Green function, G da (uj) = ((a da .; a da )), 
weighted with the Fermi-Dirac distribution function. Standard equations of motion methods 



1 14| |15 | can be used to formulate an accurate expression for G^iuj). Each level of iteration 
generates a new heirarchy of Green functions for which new equations of motion must be 
derived. To illustrate, the Heisenberg equations of motion for G da (u) generate two new 
Green functions, ((n d _ a a a ; a da )) and ((S • SdGw; a d(T )) Equations of motion for this set of 
Green functions as well as for the new Green functions that appear at this level were derived 
and solved self-consistently for the impurity density of states by invoking the Hartree-Fock 
closure. The density of states obtained at the 3rd level of iteration is sufficient to describe the 
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Kondo effect. As we have already described the T = phase, we focus on the experimentally- 
accessible high-temperature limit. The conductance 

_ -2e 2 A 



h k B T 



J f FD (u)(l-f FD (u))ImG d(7 (u + 0)du (17) 



was calculated using the standard Landauer formula J7J In Eq.(16), A = 

— Z)fc |H;| 2 5(a; — ejfc) To illustrate the role of spin-exchange, we report here the infinite U 
limit of Gd(uj) at the second level equations of motion. We find that the second level Green 
function 

G , u) = |(l-(^)) + |(S-S d ) + l(i-( Wd _ g ))-l(S-S d ) 
u - e s - E (1 - (rio-a)) u - e t - S (1 - (n d ^ a )) 

contains a contribution for the singlet and triplet states with differing spectral weights. In 
Eq. (18), the self energy is S = J2k I Vk\ 2 (^— Qc) -1 ~ — *A. This expression clearly illustrates 
that the singlet (first term) and the triplet (second term) spectral weights differ. We expect 
that the conductance into the singlet and triplet levels should reflect the asymmetry in 
the spectral weights. The conductance calculated from Eq. (17) (with the 3rd level Green 
function) is shown in Figure 1 as a function of the chemical potential. Illustrated clearly 
is the asymmetry in the conductance peaks centered at e s (the first peak) and at e t (2nd 
peak). In figure 1, Jd = —.1U. Hence, the singlet is the ground state. In the absence 
of Jdi the singlet and triplet peaks would coalesce into a single peak as in the standard 
Anderson model. The higher peaks in the conductance occur at energies 2e^ — e t + U and 
2e<i — e s +c/, respectively. The upper peaks appear inverted because for J d < 0, e s < e t . 
In the ferromagnetic regime, the triplet peak dominates and it is the neighbouring singlet 
states that lead to the asymmetry in the peak heights in the conductance as illustrated in 
Figure 2. We conclude then that spin exchange in zero magnetic field leads to peak height 
alternation in the conductance that is identical in form to the experimental || trends seen 
in the presence of a magnetic field. Ultimately, the ferromagnetic and antiferromagnetic 
cases can be distinguished by a low temperature study of the Kondo phase. 
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Figure Captions 

Figure 1: Conductance (measured in units of as a function of the chemical potential 
(measured in units of U) as computed from Eq. (17) using the third-level decoupling of the 
equations of motion for the Green function for Jd = -0.lt/ and A = 0.001c/. Each set of 
two peaks corresponds to a singlet and triplet pair, the singlet being lower in energy in the 
antiferromagnetic case. 

Figure 2: Same as Figure 1 but for the ferromagnetic coupling, Jd = 0.1U. Each set of 
two peaks corresponds to a singlet and triplet pair, the triplet being lower in energy in the 
ferromagnetic case. 
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